Shear-dependent apparent slip on hydrophobic surfaces: The Mattress Model 
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Recent experiments (Zhu & Granick (2001) Phys. Rev. Lett. 87 096105) have measured a large 
shear dependent fluid slip at partially wetting fluid-solid surfaces. We present a simple model for 
such slip, motivated by the recent observations of nanobubbles on hydrophobic surfaces. The model 
considers the dynamic response of bubbles to change in hydrodynamic pressure due to the oscillation 
of a solid surface. Both the compression and diffusion of gas in the bubbles decrease the force on 
the oscillating surface by a "leaking mattress" effect, thereby creating an apparent shear-dependent 
slip. With bubbles similar to those observed by atomic force microscopy to date, the model is found 
to lead to force decreases consistent with the experimental measurements of Zhu & Granick. 



I. INTRODUCTION 

The validity of the no-slip boundary condition is at 
the center of our current understanding of fluid mechan- 
ics. It remains however an assumption whose microscopic 
validity has been widely debated 0] . The widespread ac- 
ceptance of the no-slip condition is based on a historical 
record of outstanding agreement between theories and ex- 
periments. It is commonly agreed that the no-slip condi- 
tion results from inevitable microscopic roughness, which 
causes enough viscous dissipation to effectively bring the 
fluid to rest near the surface |^ Q . Remarkably, this 
explanation is independent of the nature of the solid and 
the liquid, contrary to ideas flrst proposed by Girard (see 

Q). 

The development of small devices has recently 
prompted a reexamination of fluid slip on length scales 
of nanometers and microns, both experimentally 1 3. Iq , 
EJlJSjJ, 10, ii, 12) and theoretically [H ElEIlll 
hj,hj,hgl. The de gree of slip is usually quantified by a 
slip length A = C/s/7, where Us is the slip velocity and 7 
is the liquid strain rate evaluated at the surface; equiv- 
alently, A is the distance below the solid surface where 
the velocity extrapolates linearly to zero f23|. In exper- 
iments, slip is usually found when the liquid partially 
wets the solid surface; measured slip lengths span four 
orders of magnitude, from molecular sizes to microns, 
and are usually shear-dependent in squeeze flow experi- 
ments, with A an increasing function of 7. In particular, 
Zhu & Granick 8] reported squeeze flow experiments, in 
which two crossed cylinders oscillate about a fixed aver- 
age distance. By measuring the viscous resistance, Zhu 
& Granick extracted the slip length over a wide range 
of oscillation amplitudes. These experiments lead to the 
largest shear-dependent slip lengths yet (up to ~ 2 //m). 

The origin of this large shear dependent slip is hereto- 
fore mysterious. Nanobubbles have recently been ob- 
served on hydrophobic surfaces, using atomic force mi- 
croscopy, with typical thickness h ^ 10 nm, typical 
radius i? ~ 50 — 100 nm and high surface coverage 
HI m m 113. Ahhough the origin of these bub- 
bles is unclear and skepticism remains in the community 
about their existence, they have been often invoked as 



apossible origin o f the so-called hydrophobic attraction 
|25ll26ll27ll23 . l29l | and their existence points to a possible 
picture for such large slip [T^l2^ . 

It is well known that there is in general a non-zero ve- 
locity at a liquid-gas interface, and therefore it is natural 
to wonder whether the existence of such a gas layer at 
the solid surface is sufficient to explain the experiments. 
When a fluid of viscosity r]i adjoins a layer of fluid of 
thickness h with smaller viscosity 772, the discontinuous 
strain rate at the fluid-fluid interface results in an appar- 
ent slip with slip length 



A = M^-1 



(1) 



Choosing 771/772 = 50 appropriate for a gas- water inter- 
face leads to slip lengths as large as 500 nm. This esti- 
mate is however independent of the interfacial shear and 
therefore unable to explain the squeeze flow experiments; 
it also overestimates the slip length in the case of bubbles, 
as is discussed in section imi Note however that similar 
arguments are consistent with data from pressure-driven 
flow experiments where reported slip lengths to date are 
essentially shear-independent [lol 111 [ial30l Isil Is^ . 

In this article, we will assume bubbles exist on hy- 
drophobic surfaces and will calculate their dynamic re- 
sponse to an imposed oscillatory shear. In an oscillatory 
squeeze flow experiment fl, we argue that the pressure 
fluctuations in the fluid cause the bubbles to act as a 
"leaking mattress", with both compression and dilation 
of the gas in the bubble, as well as diffusion of gas into 
(and out of) the bubble. As the solid sphere oscillates, 
this periodic in-phase response of the bubbles sizes re- 
duces the amount of liquid necessary to be squeezed out 
of the gap and thereby the force on the moving sphere, 
creating an apparent slip. Our calculations indicate that 
the magnitude of this apparent slip is consistent with the 
observations of Zhu & Granick. We present the details of 
our model in the next section and discuss the comparison 
with the experiment in section ITTll 
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II. INFLUENCE OF BUBBLES ON FORCE 
MEASUREMENTS 

A typical oscillatory squeeze flow experiment is shown 
in Figure n A sphere of radius a oscillates with velocity 
Vs in a viscous liquid of viscosity 77 at a distance D of 
a planar surface (equivalently, the surfaces can be two 
crossed cylinders). The two surfaces are assumed to have 
the same physico-chemical properties. If no bubbles are 
present and the no slip boundary condition is satisfied on 
both surfaces, the lubrication force opposing the motion 
of the sphere is given by the Reynolds equation 

Fit)=e,-F{t) = -^^Vs^F,^^. (2) 

If however flow occurs on the surfaces with a slip length 
A, the viscous force is decreased by an amount fsUp |33| 
given by 
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(3) 



Equation @ is used experimentally to infer effective slip 
lengths: the experimental viscous force i^exp is compared 
to the theoretical no-slip result i^iub and any difference is 
interpreted as fluid slip, with a slip length A correspond- 
ing to /slip = i^oxp/^lub- 

Let us now assume that the solid surface is covered 
with a percentage (j) of identical gas bubbles (Figure^, 
and determine how the bubbles modify the dynamic re- 
sponse. Although this assumption has been made by pre- 
vious authors 12^ . |2^ I27I I2M |29| , the physical mecha- 
nism responsible for such bubbles is unknown. Simple es- 
timates indicate that small bubbles are short lived in so- 
lution [s^. However, stable bubbles could arise from any 
number of possibilities that are known to prolong bubble 
lifetimes, including surfactants, surface heterogeneities, 
or local supersaturation of dissolved gas In this 

paper, we are interested in understanding whether the 
dynamic response of hypothesized bubbles is sufficient to 
rationalize slip experiments. 



A. Total force 

The total force F{t) resisting the oscillatory motion of 
the area S of the sphere has two components: (1) a vis- 
cous lubrication force Fh , due to hydrodynamic pressure 
fluctuations and acting on an area (1 — (j))S and (2) an 
elastic bubble force F},, due to pressure fluctuations in- 
side the bubbles and acting on an area (j)S. The total 
force is therefore given by 



where 



F{t) = (1 - cl>)Fh + <l>Fb, (4) 



Fh^ ip- Po)S, Fb = {pb - Peq)S. (5) 




FIG. 1: Typical squeeze flow experiment: a solid sphere of 
radius a is oscillated in a liquid at a distance D <C a of a 
smooth solid surface with amplitude d <^ D and frequency 
uo. The surfaces are covered by microscopic gas bubbles of 
contact angle 9 and radius of curvature Rq. The set of bubbles 
is approximated by a gas layer of time evolving thickness h{t). 



Here p and pb {po and Peq) denote the (equilibrium) pres- 
sures in the liquid and the bubbles respectively. More- 
over, since D <^ a, the surface S is given hy S ~ naD. 



B. Lubrication force 

Let us first calculate the hydrodynamic force Fh- 
The presence of bubbles modifies equation Q in two 
ways. First, flow occurs over a distribution of bubbles 
located on an otherwise no-slip surface, so the viscous 
force is reduced by an amount fsUp given by equation 
(PI, where A is the appropriate effective slip length for 
flow over a distribution of bubbles 0, [s^ . 

Second, the size of the bubbles changes in time in re- 
sponse to pressure fluctuations in the liquid. This volume 
effect will modify the amount of liquid necessary to be 
squeezed out of the gap at each cycle of the oscillations, 
hence the viscous force. Consequently, bubble dynamics 
has to be subtracted from the forcing Vs and the hydro- 
dynamic lubrication Fh force is now given by a modified 
Reynolds equation 



Fh = 



slip " 



15~ 



Vs 



dh\ 



(6) 



where h is an average bubble thickness on each surface 
and the factor 2 accounts for the fact that each surface 
is covered with bubbles. 



C. Rate of change of bubble height 

In order to calculate dh/dt in equation ©, let us now 
consider the dynamics of the bubbles. We assume the 
bubbles are undeformed by viscous stresses and remain 
spherical, with radius of curvature R(t) and interior angle 
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(tt — 9) (see Figure^. We neglect interactions between 
bubbles. We expect h{t) to depend explicitly on the forc- 
ing on the bubbles, i.e. Ff,. 

At the small frequencies typical of squeeze flow exper- 
iments (1-100 Hz), the gas is isothermal, so the pressure 
in the bubble changes via the ideal gas law 



Pb{t)V{t) Pe,VQ 



m{t) 



mo 



(7) 



where V and m denote the volume and the mass of a 
single bubble. The average thickness of the gas layer is 
defined as h(t) = nV{t), where n is the number density 
of bubbles on the solid surface, so that equation (|7J) can 
be rewritten as 



Pb{t)h{t) _ Peghp 
m{t) niQ 



(8) 



Combining the time-derivative of © with Fi, = {pi, — 
Peq)S and linearizing around {pb,h,m} ~ {peq,hQ,mQ\, 
we obtain the equation for the rate of change of the mean 
bubble height h 



dh ho dm ho dFh 
dt mo dt PeqS dt 



(9) 



We thus have that the rate of change of h is the sum of 
a rate of change governed by gas diffusion plus a second 
contribution due to the gas compressibility. 

We now consider the rate of gas diffusion from the 
bubble. In our model for the oscillatory squeeze ffow 
experiments _2., ,6, 8j , bubbles lose mass by both vertical 
diffusion across the liquid layer and radial diffusion along 
the apparatus; because of the scale separation Z3 <C a, 
these two processes require separate treatment. 

Let us first consider the case of vertical diffusion. Since 
for most common gases, k ^ 10^^ m^/s, the vertical 
Peclet number Pcy = D^oj / k is much smaller than unity: 
on the experimental time scale w^"'^, the bubbles are ap- 
proximately in instantaneous vertical diffusive equilib- 
rium. The dissolved gas concentration above the bub- 
ble is therefore uniform throughout the liquid gap and 
is given, in the case of small amplitude oscillations, by 
Henry's law c = PbCao/Peq- The mass m of gas necessary 
to fill the liquid gap at this concentration is equal to the 
gap thickness {D — 2h) times the change in concentra- 
tion Coo(j>b/Peq — 1) timcs the area in the liquid which is 
influenced by the bubble, i.e. 1/n . Linearizing around 
h ^ ho and combining with Peq/cao — Po/^o, we get that 
m is proportional to the bubble force Ft 



co{D - 2ho) 
npoS 



Fh- 



(10) 



Consequently, the total rate of change of m is given by 
dm dm,. co{D — 2/io) dFf, 



dt 



dt 



npoS 



dt 



(11) 



where drur/dt is the rate of change in the bubble mass 
governed by gas diffusion in the (slowly varying) radial 



direction of the apparatus; let us now evaluate this con- 
tribution. 

In contrast to the vertical case, the radial oscillatory 
Peclet number, Per = L^lo/ n = qDuj/k, is of order unity 
or larger, so that radial diffusion has to be accounted 
for explicitly. Assuming the dissolved gas is in vertical 
diffusive equilibrium, the time rate of change of the mass 
of a gas bubble drrir/dt is given by a ffux integral on the 
bubble surface Sb 



drrir 
"dT 



(12) 



where the assumption od spherical cap bubble implies 
that I{e) = 7r(2(7r -6)+ sin 29) /2. 

In general, the radial concentration of dissolved gas 
c(r, t) verifies an advection-diffusion equation with shear 
dependent diffusivity [s^l- However, for the small ampli- 
tude oscillation in j^] , both advection and Taylor disper- 
sion are neghgible, and c{r,t) satisfies a pure diffusion 
equation. We finally approximate radial concentration 
gradient by a simple linear law dc/dr k, (cqo — c)/Lr 
where ~ {kIuj)^I'^ w L/PeV^ is the typical (shear 
dependent) radial gradient length scale. Equation l|12(l 
together with Henry's law leads therefore to a linear re- 
lation between the rate of change rrir and the bubble force 



dm,. kRII{9)cq 



dt 



PoSLr 



Fb- 



(13) 



Combining mo = poho/n with l^, 111|) and 113|) . we 
finally obtain that the mean bubble height h satisfies the 
differential equation 



_ = .k,Fb-k, — , 



where (fci, ^2) are given by 
nKRoI{0)co 



PoPoSLr 



coho f Coo{D-2ho) 

K2 = I 1 



CooPO'5' 



(14) 

(15) 
(16) 



D. Bubble force 

We finally need to calculate the bubble force Fb in or- 
der to close the system of equations and (|14|) . This 
in general requires understanding the (equilibrium or 
nonequilibrium) mechanism responsible for the presence 
of these long-lived bubbles. However, we can bypass this 
unknown physics by assuming without loss of generality 
that the pressure fluctuations in the bubbles and that in 
the liquid are proportional 



Apb = a Ap, 



(17) 
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where a is an unknown constant. Hence, bubble and 
hydrodynamic forces are proportional Fj, — aF^ , and the 
total force on the sphere can be expressed as 



F 



ad}) Fh = 



1 — + acf) 



Fb. 



(18) 



E. Total force on the sphere 

We can now combine H14() and (|18|l to express the 
total force opposing the motion of the sphere F. We 
obtain 



F{t) = --{l~-c^ + a^ 



Vs + 2k,Fb 



2k/-^ ] , (19) 



where(5 12Trria^ fsUp / D . Since F and Fi, are related by 
equation (|18|l . (|19|l can be transformed into an ordinary 
differential equation for F, 

F{t)^-^(vs + 2akiF + 2ak2^y, (20) 

where Vs = (1 — + a(f))Vs- For an oscillating sphere 
velocity Vs = d (dsmuit) /dt, the periodic solution to 
(pHl is given by 

Fjt) ^ (l-0 + a0)(l + fofci) 
Fi,b ^'"P(l + fofci)2 + (J^afc2)2 ^ ^ 

SLuak2 

X I costjiH — —sm ujt 

1 + oafci 



F. What is the value of a? 

By comparing our model pifl to the results of squeeze 
flow experiments, we find that the only choice consistent 
with available data at large separation distances is a w 1. 

To see this, consider equation H21|l in the limit of large 
separations between the sphere and the planar surface 
D. Since equation (jS)) shows that fsiip ~ 1 when D is 
large, we get 6 ~ D~^. Moreover, S ^ D so that, from 
((T^ and we obtain Ski ~ F)~^ and 6k2 ~ D~^. 

Consequently, in the limit of large separations, we obtain 
that the ratio of the measured force, out-of-phase with 
the sphere displacement, to the expected no-slip force Q 
is given by 



lim 

D^oo 



F \ _ ^ 

-Flub/ 



(22) 



Within experimental errors, this ratio is always measured 
to be unity 0, Q, IE IS Q , *■ e- the expected lubrication 
no-slip force is recovered for large separation distances. 
We therefore need 1 ~ (j> + acf) « 1 or a « 1. 

We emphasize that this conclusion is reached because 
we assume that the model presented in ^ is the major 
physical mechanism responsible for the force decrease ob- 
served in experiments such as 13- 



G. Final formula for the force ratio 

We obtain from (|21|1 that the ratio /* of the peak force 
out-of-phase with the sphere displacement to that ex- 
pected with no-slip and no bubbles (|2Jl is given by 



/slip 



1 



5ki 



{uj5k2Y_ 

l + 5ki 



(23) 



The "leaking mattress" model therefore leads to an ap- 
parent slip effect, of dynamic origin. The effect is shear- 
dependent through the frequency dependence in H23|l : 
higher frequency and therefore higher shear shear rates 
lead to a larger apparent slip, in agreement with 0.SI3- 
The model was derived under the assumption of small 
amplitude oscillations d, which consequently does not ap- 
pear in the final formula for /*. 

Note that in the limit of high frequencies the force ratio 



(|^ becomes /*(w 
cies J*{u} — > 0) 



oo 



whereas at low frequen- 



jqjj^. In the low frequency limit /* 
is independent of frequency and depends only on sepa- 
ration distance. Furthermore, equation H23I) implies that 
the apparent slip effect increases with the fluid viscosity, 
in agreement with experiments 0, ^3 . It also increases 
with the size of the sphere a, which might account for the 
large slip lengths reported in |3] (cm-size spheres) as op- 
posed to other squeeze flow experiments (usually /im-size 
spheres). The model predicts that the measured overall 
apparent slip length is therefore not only a solid/liquid 
property but depends on the system size • 

We finally note from (|21|l that the total pressure force 
on the sphere F{t) also contains a non-zero component in- 
phase with the displacement of the sphere, and therefore 
out-of-phase with its velocity. If we denote by g* the 
ratio of this in-phase response to the expected no-slip 
no-bubbles out-of-phase response, we obtain 
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6u!k2 



{I + skiy + iSijk2y 



r /slip 



Sujk2 
l + Ski 



r 



(24) 



Equation 1)24(1 is a prediction of the effective elasticity 
provided by the bubbles to the surface, which would oc- 
cur in addition to other in-phase contribution such as in- 
termolecular forces, and is experimentally testable. The 
values of the in-phase responses of the forces were unfor- 
tunately not reported by Zhu & Granick 



III. COMPARISON WITH EXPERIMENTS 

We present in this section a quantitative comparison 
of our model with the experimental results of Zhu & 
Granick in the case of deionized water, namely the four 
sets of data presented in Figure 2 of ref • The macro- 
scopic water/solid contact angle in this case was 110° the 
sphere radius was a = 2 cm and we assumed that the liq- 
uid was saturated with O2 at 25°C and 1 atm (cqo = cq), 
for which po = 1.28 kg/m'^, cq = 8.3 x 10"'^ kg/m^ and 
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FIG. 2: Comparison between the experimental data of Zhu & 
Granick (2001) and the dynamic model 12311 with iio = 10 nm 
and = 99%. (a): Small amplitude experimental data; (0): 
measurements for d — 0.5 nm, ui — 1 Hz; dashed-dotted line: 
model for 9 = 177°; (l>): measurements for d = 1.6 nm, uj = 
10 Hz; solid line: model for 9 = 132°. (b): Large amplitude 
experimental data. (♦): measurements for d = 6 nm, u — 1 
Hz; dotted line: model for 9 — 168°; (►): measurements for 
d = 6 nm, u) = 10 Hz; dashed line: model for 9 = 90°. 



K — 2 X 10^^ m^/s. As a matter of comparison, we have 
also summarized in Table ^ the experimental results of 
[2lll2^[23ll2^ on the typical size, distribution and mor- 
phology of bubbles observed by atomic force microscopy. 

The "leaking mattress" model we have presented in 
the previous sections has three free parameters which we 
fit to the experimental data: (a) the area fraction of the 
bubbles on the surface, < ^ < 1, (b) the size of the 
spherical bubbles, described by their radius of curvature 
Rq and (c) the microscopic contact angle 6 at the bub- 
ble level, which significantly differs from the macroscopic 
contact angle because of both intermolecular forces at 
the nanometer scale. Note that <j) is related to the area 
fraction n by the formula (jy—nnRQ sin^ 9. 

Furthermore, in order to present a meaningful fit to 
available data, we require that in each experiment the two 
layers of bubbles fit in the gap between the sphere and the 
plane for all separation distance. This is a geometrical 
constraint written as 2ho — 2Ro{l + cos9) < min(£') |Q 

The model H23|) can be well fit to the experiments Q 
with appropriate parameter choices. The best fits are 
obtained when we choose Rq « 10 nm. This is illustrated 
in Figure [21 where the fits are compared the force ratio 
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FIG. 3: Comparison between the experiment from [g with 
d = 1.6 nm and = 10 Hz (l>) and the model for different 
bubble sizes, contact angles and surface coverage, (a): Influ- 
ence of bubble size; model with surface coverage = 0.99, 
contact angle 9 = 150° and bubble sizes _Ro =1 nm (dotted 
line), 25 nm (solid line) and 50 nm (dashed line), (b): Influ- 
ence of contact angle; model with surface coverage (/> = 0.99, 
bubble size J?o = 10 nm and contact angles 9 — 120° (dashed 
line), 140° (solid line), 160° (dotted line) and 170° (dashed- 
dotted line), (c): Influence of surface coverage; model with 
bubble size Rq = 10, contact angle 9 — 120° and surface cov- 
erage = 0.1 (dashed-dotted line), (p = 0.25 (dotted line), 
(l> = 0.5 (solid line) and (j) — 0.99 (dashed line). 



from the model to the small (a) and large (b) amplitude 
data from Zhu & Granick j^; the values of the angles 9 
were chosen for each curve to be the best in a least-square 
sense and <j) — 99%. As expected from the linearity of our 
model, the fit to the low amplitude data of is better 
than that obtained for oscillations of larger amplitude. 

We explore the infiuence of the three parameters of 
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our model {(j), Rq, 6) in Figure|3t-c for the measurements 
from with d — 1.6 nm and u — 10 Hz. 

We first find that the results of our model depend 
weakly on the bubble sizes: the results on Figure ISt 
are consistent with the experimental data for a large 
range of bubble sizes, i?o ~ 1 to 50 nm. These sizes 
are agreement with the experimental evidence of bubbles 
in [21^, 22, 2dj 2A] as summarized in Table HI although 
somewhat smaller. As a matter of comparison, the data 
in show large standard deviation (up to 70%) for the 
bubble area. 

As a difference, we find that the results of our model 
depend on both the assumed microscopic contact angle 9 
and coverage of the surface by the bubbles (f>. We observe 
variations in the contact angles leading to best fit to the 
four experiments (Figure EJ and also note that we obtain 
a departure from the best fit when the angle is chosen 
to be significantly different (Figure I^Js). In three experi- 
ments out of four, we find that the microscopic contact is 
larger than the macroscopic contact angle 110° character- 
iziiig the wetting of deionized water on the surfaces used 
in j^]. This result is consistent with the data in Table HI 
where, in all cases, bubbles were found experimentally to 
be fiat with microscopic contact angles larger than the 
macroscopic wetting angles. The fourth set of data from 
is found to be consistent with a microscopic angle of 
about 90°. Although this is different from the data in 
pH 123 . 123L 01 1 it is consistent with theoretical studies 
which show that intermolecular forces lead to microscopic 
contact angles which are always closer to 90° than their 
macroscopic counterpart [s^. Furthermore, we note that 
electrical effects are known to have significant impact on 
contact angles of bubbles and drops 39, 40]. 

Finally, we find that our model is consistent with the 
experimental data when the surface coverage is assumed 
to be large and almost equal to 100% (see Figure Et) . 
This result compare well with the available data on bub- 
bles where, in three out of four studies [l^ 113,13' the 
bubbles were found to cover almost entirely the solid sur- 
face. As a difference, the pictures in show bubbles 
with lower surface coverage. We also note that our pre- 
vious study of slip in pressure-driven flow experiments 
lead to a similar conclusion: in order for surface-attached 
bubbles to be responsible for the measured effective slip, 
surface coverage of almost 100% was necessary [T^ . 



IV. CONCLUSION 

We have explored in this paper the consequences of the 
presence of nanobubbles on the surfaces where squeeze 
flow experiments are performed. We have shown that, 
within the framework of a simple stabilizing model, the 
time-dynamics of bubbles always leads naturally to a 
shear-dependent decrease in the measured viscous force 
by a "leaking mattress" effect. The effect was found to 
increase with viscosity of the fluid and the size of the 
sphere, in agreement with earlier experimental results. 



We emphasize that this mechanism is of dynamic ori- 
gin, and is not a consequence of the microscopic slip at 
the bubble surfaces; in particular, we argue that this is 
why shear-dependent slip length have not been reported 
by investigations of slip in pressure-driven flow experi- 
ments to date, where no oscillatory pressure is present to 
trigger an effect similar to the one proposed here. Also, 
the mechanism we propose should also apply to squeeze 
flow experiments performed on super-hydrophobic sur- 
faces such those reported in jl^] with small air bubbles 
trapped on fractal surfaces (see also 0)- 

Assuming the presence of bubbles, the calculations on 
the model have been performed with several simplify- 
ing assumptions and, in particular, additional contribu- 
tions to the sets of coefficients (fci, fc2) could come from 
bridging bubbles, large amplitude oscillations of the solid 
sphere or bubble interactions, deformation or displace- 
ment on the solid surface. 

We have then presented a comparison between the re- 
sults of our model when applied to the experiments of 
Zhu & Granick. We found that our model gives results 
which are in agreement with the force decrease measured 
experimentally, for bubble features which are consistent 
with available experimental data on nanobubbles (bubble 
size Rq lOnm, large microscopic contact angles, large 
surface coverage). Finally, a formula has been proposed 
for the (additional) effective elasticity provided by the 
bubbles to the solid surface. 

We note that our study does not rule out the possi- 
bility of bubbles with dynamically selected sizes. It has 
been reported experimentally in [4^ that the jump-in 
distance between two hydrophobic surfaces in water, be- 
lieved to be due to the presence of bubbles, depended 
on the history of the sample; performing the experiment 
several times lead to changes in the jump-in distances 
over time which was found to remain constant only after 
a few periods. A similar scenario could be envisioned in 
experiments such as 

To conclude, we present a simple prediction based on 
the results of our model. If a squeeze flow experiment was 
performed with two different surfaces, say a hydropho- 
bic plane and a hydrophilic sphere, force ratio measure- 
ments displaying shear-dependant results should be able 
to test whether the ideas put forward in this paper are 
valid. Indeed, if the force decreases was really due, not 
to bubbles, but to a change in the hydrodynamic bound- 
ary condition for flow past the hydrophobic surface, the 
maximum force decrease one could expect to obtain is 
1/4 for the case of a perfectly slipping surface (see [s^ 
for the calculation; this result can also be found by sym- 
metry about the plane where slip occurs) . If alternatively 
the measurements are due to a "leaking mattress" effect 
similar to the one we propose here, equation ^ should 
also apply (with different prefactors) and therefore so is 
equation (|23|l : consequently, force ratio smaller than 1/4 
should be obtained in this case. This proposition for an 
experiment, together with the prediction for the in-phase 
response of the force (|24|) . would allow our model to be 
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tested experimentally. 
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Ishida et al. [21J Tyrrell & Attard [22] Tyrrell & Attard [23] Steitz et al. [24] 

Projected area (nm^) 3.3 x 10^ 4 - 6 x 10^ 4 - 7 x 10^ 2 - 11 x 10^ 

Height ho (nm) 40 20 - 30 20 - 30 < 18 

Radius of curvature i?o (nm) 1300 ~ 50 40 - 60 30 - 100 

Surface coverage <j> ~ 20% ~ 100% ~ 100% 89% 

Macroscopic contact angle 110° 101° 101° > 90° 

Microscopic contact angle 166° ~120° 117° -130° 130° -147° 



TABLE I: Summary of experimental data on nanobubbles as 
found in 21, 22, 2^ |2^ by atomic force microscopy: pro- 
jected area of each bubble on the solid surfaces, height above 
the surfaces ho, radius of curvature Ro, surface coverage (j), 
macroscopic and microscopic contact angle. The radius of 
curvature and microscopic contact angles were inferred from 
the other data assumming spherical cap nanobubbles. 



